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ABSTRACT 

I discuss a new gravitational potential, <I>(r) oc (rj, 1 + r n ) — 1//n , for modeling the mass distribution 
of spherical systems. This potential has a finite mass and generates a density profile with inner slope 
2 — n. A gas embedded in this potential has hydrostatic temperature and gas density distributions that 
are elementary functions of n, greatly simplifying the task of measuring the slope from X-ray data. I 
show that this model is successful in describing the rising temperature profile and steep gas density 
profile often seen in cooling flow clusters. An application to the Abell 478 cluster of galaxies yields an 
inner slope 2 — n = 1.0 ± 0.2 (90%), consistent with the inner regions of collisionlcss dark matter halos 
first simulated by Navarro, Frenk, and White. The potential is also useful for cluster dynamics: it is a 
generalization of the familiar Hernquist and Plummer potentials, and because it is invertible, it allows 
for easy analytic calculation of particle phase space distribution functions in terms of n. 
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1. INTRODUCTION 

The X-ray emitting medium in clusters of galaxies is 
gravitationally heated, and therefore a tracer of dark mat- 
ter. The current generation of X-ray satellites finally pro- 
vides the spatial and spectral resolution necessary to con- 
strain not only the total amount but also the distribution 
of matter in clusters. The shape of the mass distribution 
is important, because it can be usefully compared with 
N-body simulations of structure formation, allowing con- 
straints on dark matter physics. For example, the inner 
slope of the dark matter density, 
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can distinguish between collisionlcss (y > 1) and self- 
interacting [y « 0) varieties of dark matter (Nakano & 
Makino 1999; Moore et al. 1999; Burkert 2000). 

Direct measurements of v from X-ray data are rare, how- 
ever, for two reasons. The first is that most mass models 
in common use — e.g., the Navarro, Frenk, & White (1997; 
NFW) profile p oc (ro+r) 2 jr or the nonsingular isothermal 
sphere — do not have an adjustable inner slope v, and even 
if they do, the gas and temperature distributions they gen- 
erate are not expressible in terms of elementary functions 
of r and v (Markevitch & Vikhlinin 1997). This inconve- 
nience makes fitting for v an expensive and difficult task. 

More importantly, our understanding gas physics in the 
cores of the galaxy clusters closest to equilibrium leaves 
something to be desired. In the ROSAT era, such clus- 
ters were observed to have cool cores and cooling times 
much shorter than the Hubble time; for these reasons they 
were thought to possess a quasihydrostatic cooling flow 
(Fabian 1994; Peres et al. 1998). However, the Chandra 
and XMM satellites have in general not detected the mul- 
tiphase plasma required by cooling flow models (Bohringer 
et al. 2001; Johnstone et al. 2002; Lewis et al. 2002). It 
is therefore important at least to ask whether true hy- 
drostatic models are capable of describing the properties 
of the gas — specifically the rising temperature profile — 
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outside the central « 10 kpc, where the gas is often dis- 
rupted by a central radio source (Blanton et al. 2001). 

In this Letter I consider a unified solution to the above 
two problems. First, I discuss a mass distribution with 
adjustable inner slope which does yield hydrostatic gas 
and temperature profiles in terms of elementary functions. 
Then, I show that together with an exponential pressure 
distribution, these hydrostatic models can reproduce the 
rising temperature profile and the steep density profile in 
the rich cluster Abell 478. 

2. THE POTENTIAL 



Consider a potential of the form 

GM tot 



$(r) 
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Where M tot is the total mass of the system. The cumu- 
lative mass distribution and density profile that generates 
this potential are 

/ n \ (l+n)/n 

M ( r ) = M tot -zr—Z (3) 
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The requirement that M(0) — yields the lower bound 
n > 0, while the requirement that the density never in- 
crease with radius leads to an upper bound n < 2. This 
mass profile therefore corresponds to a broken power law 
with inner slope v = 2 — n and outer slope 3 + n. The 
n = 1 case corresponds to the Hernquist (1990) model, 
n = 2 is a Plummer sphere, and n — 1/2 yields a profile 
similar to that of collisionless dark matter halos in recent 
N-body simulations (Ghigna et al. 2000). 

The potential is invertible, and thus the density can be 
expressed as a function of the potential. If we define a 
dimensionless potential variable ip = —^tq/GM then 

M tot (n + l) V™ +3 
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This form of the density can be used to compute phase 
space distribution functions using the methods outlined by, 
e.g., Cuddeford (1991), and thus facilitate measurement of 
the mass profile through equilibrium particle dynamics. 

3. X-RAY ATMOSPHERES OF DARK MATTER HALOS 

Here I investigate the luminosity and temperature dis- 
tribution of a hot, X-ray emitting plasma embedded in 
the mass distribution given by (3) and governed by the 
equation of hydrostatic equilibrium: 

VP = -p 9 V$, (6) 

where P is the gas pressure and p g is the gas density, which 
in general is different from the total matter density traced 
by <E>. With the assumption that the plasma is spherically 
symmetric, the equation becomes 



1 dTp g 
Pa dr 



dr 
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where T = P/p g is the thermal energy per unit mass. For 
an ideal gas in a cluster environment, T = kT/p,m p , where 
T is the temperature, m p is the proton mass, p = 0.61 is 
the mean molecular weight, and k is Boltzmann's constant. 

3.1. Polytropes and the f3 -model 

Consider the possibility that the gas is a polytrope, 
P oc p 1 , where 7 is the polytropic index. In this case, 
the potential yields 
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3(7 ~ 1) r T 
An important property of this model is that the n = 2 
potential reproduces the familiar /3- model, the gas density 
profile most widely used by X-ray astronomers. While 
most literature assumes that the /3-model describes an 
isothermal gas, in the present context it corresponds to 
the temperature profile in equation (8). 

Furthermore, the model suggests an important connec- 
tion between the polytopic index 7 and the parameter. 
It is clear from equation (9) that, for the gas mass to 
be finite, we must have > 2/3. On the other hand, for 
< 2/3, the polytropic index 7 > 5/3, and the plasma be- 
comes subject to convective instabilities (Nevalainen et al. 
1999). Thus, two important physical criteria for the valid- 
ity of the polytopic model above — convergence of the gas 
mass and lack of convective instabilities — are fulfilled by 
the very same constraints on 0. 

3.2. A Hydrostatic Alternative to Cooling Flows 

The polytropic models have two marked disadvantages: 
they require that the X-ray emitting medium have the 
same characteristic radius, r , as the total matter distri- 
bution, which need not in general be true; they also do not 
reproduce the temperature structure of the clusters clos- 
est to equilibrium. In many of the most regular clusters 
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of galaxies, the temperature is observed to rise with dis- 
tance from the center, leveling out or falling after a peak. 
However, it is clear from equation (8) that the tempera- 
ture in the polytropic model decreases monotonically with 
radius, and is therefore unable to reproduce the rise in 
temperature observed in many clusters. 

For this reason, I introduce a new hydrostatic pressure 
model that reproduces many of the properties observed 
in clusters with rising temperature profiles. Consider an 
exponential gas pressure profile: 

P = P e- {r ' ra ^l b (11) 

where Pq is the central pressure, and b 1/>a ro is the scale 
height of the pressure distribution. Together with the 
equation of hydrostatic equilibrium and equation (2), this 
yields the following gas density and temperature profiles: 
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In order for the gas density to decrease with the radius, we 
must have n > a. In this case, the temperature profile be- 
haves as an increasing power law of index n — a, achieving 
a maximum at r max = r [(n — a)/(l + a)] 1 /" before falling 
to zero at infinity. When n = a, the temperature peak 
is at r = 0, giving us again a monotonically decreasing 
temperature profile. 

In Figure 1, I show that this model fits the temperature 
and the emissivity profile of Abell 478 as observed by the 
Chandra satellite. 2 The deprojected data are from Sun 
et al. (2002); all analysis and model fitting is mine. The 
minimized six-dimensional merit function is 
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where N n is the total number of density data points rij, 
and Nt is the total number of temperature data points 
Tj. The minimized value is 39.0 for 49 degrees of freedom, 
indicating that the fit is of a high quality. 

There are 6 total parameters, and therefore 6 x 5/2 = 15 
unique pairs of parameters. I derive marginalized joint 
probability distributions for all these pairs via a four- 
dimensional integration of the Bayesian probability dis- 
tribution exp(— x 2 /2), with uniform prior. Six of these 
marginalized distributions are shown in Figure 2. To ob- 
tain closed confidence contours, it was necessary to substi- 
tute, without loss of generality, the parameters r e = 6 1//a ro 
and M0.5 = M(< 0.5 Mpc). These substitutions are sensi- 
ble, as r e is just the exponential scale height of the gas, and 
M5Q0 is the mass within a spherical radius of 500 kpc — 
about the extent of the Chandra field of view. 

Finally, in order for the previous discussion to be phys- 
ically self-consistent, three criteria must apply. First, the 
gas density must never exceed the total matter density. 
This corresponds to the requirement that the baryon frac- 
tion pg/p < 1 everywhere throughout the cluster. Sec- 
ond, the gas must not experience convective instabilities 
(-G?lnp g /<imr < 5/3). Third, the inner slope of the dark 
matter density pdm = P — Pg must be checked against the 
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inner slope of the total density p. Figure 1 shows that all 
these criteria are abundantly fulfilled. The average baryon 
fraction within 0.5 Mpc is 0.08 ±0.07 (at 68% confidence), 
and the radial variation of it within that range is equally 
well-constrained. The effective polytropic index is always 
< 5/3, and the slope of the dark matter density profile 
matches that of the total profile. Hence the model as ap- 
plied to Abell 478 is self-consistent. 

Of course, any hydrostatic model has to contend with 
the fact that the cooling time remains quite short. Figure 
1 shows that the cooling time on the average is smaller 
than the Hubble time in the region observed by Chandra, 
in agreement with Sun et al. (2002); the gas should radiate 
away all its energy in under ten billion years. However, the 
failure of Chandra and XMM observations to detect the 
temperature distribution predicted by cooling flow models 
(Lewis et al. 2002; Peterson ct al. 2001) has lead some to 
infer that an unknown heating mechanism must balance 
the cooling (Loeb 2002; Bruggen & Kaiser 2002). While 
the details of this balancing system are beyond the scope 
of this Letter, its net effect may be to make the hydrostatic 
solutions empirically valid in at least some cases. 



4. CONCLUSION 

I suggest a new potential, $ oc (r™ + r% )~ 1 ^ n , for mod- 
cling the matter distribution of galaxy clusters. This po- 
tential has an inner density cusp of slope 2 — n, and is 
particularly useful because (1) it is a generalization of po- 
tentials in common use; (2) it is invertible; and (3) a gas 
embedded in it can have density and temperature distribu- 
tions that are elementary functions of n. As an example, I 
successfully apply a hydrostatic model to published Chan- 
dra observations of the galaxy cluster Abell 478, which is 
generally believed to contain a cooling flow. The resulting 
inner slope 2 — n = 1.0 ± 0.2 (90% confidence interval) 
is well constrained, and consistent with N-body simula- 
tions in which massive halos are composed of collisionlcss 
matter. 

This work demonstrates that cold cores and steep den- 
sity gradients are compatible with hydrostatic models in 
all respects except for the cooling time. 
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the Smithsonian Astrophysical Observatory for and on be- 
half of NASA under contract NAS 8-39073. 
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Distance From Cluster Center (kpc) 

Fig. 1. — Fits of the present model to the temperature and electron density distribution of Abell 478. The data are taken from Sun et al. 
(2002); the model fits are mine, (a) Dcprojected temperature profile; (b) deprojected electron density. Shown also are 68% confidence limits 
on the following quantities as a function of radius: (c) the cooling time; (d) the slope of the dark matter density v = dlnp(j m /dlnr; (e) 
the effective polytropic index d In P/d In p g ; and (f) the baryon fraction (solid lines) as well as on the average baryon fraction with 0.5 Mpc 
(dashed horizontal lines). 
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- Marginalized Bayesian probability distributions for 6 of the 15 unique parameter pairs. Shown are the 68%, 95%, and 99.7% 
contours. 



